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ARTIFICIALLY SHOCKED GRAPHITE' 


THOMAS P. SCIACCA AND MICHAEL E. LIPSCIIUTZf 


Sample.: of graphite, partly converted to diamond by shock, and meteoritic graphite were 
studied by transmission election microscopy. The shocked graphite remained in the form of 
single crystals and was mildly deformed while the diamond which had formed was anhedral end 
polycrystalline. The meteoritic graphite contained minor amounts of troilite (FeS) and elemental 
sulfur but no detectable carbon phase other than graphite. 


As an adjunct to an x-ray crystallographic study 
of meteoritic diamonds (/) we have studied the 
morphology of the phases present in meteoritic 
graphite, and in graphite which has been artifi- 
cially shocked at pressures high enough (300 kbar) 
to produce diamond ( 2 ). 

Samples of the starting material, high-purity 
graphite (3), as well as the resulting shock- 
product, and graphite from the interior of a 1.1-g 
graphite nodule from the Canyon Diablo iron 
meteorite were crushed to 100 mesh. Slurries in 
ethanol were then placed in an ultrasonic gene- 
rator for 5 minutes to disperse the sample uni- 
formly. The suspensions were placed on 
Formvar, collodion, or carbon substrates for 
examination in a Hitachi HU-1 1 electron micro- 
scope. Phases were identified by electron diffrac- 
tion of selected areas with gold and palladium as 
standards. 

Figure la shows a typical graphite-diamond 
intergrowth from the artificially shocked graphite. 
All of the diamonds formed in this manner were 
anhedral and polycrystalline. The graphite 
flakes seemed mildly deformed although they 
remained single crystals. Figure 16 shows a 
typical graphite flake from the starting material. 

In both the starting material and shock product 
the only other phase that we found was in the 

‘Published in Seine*. 14S (3030): )04»-1080. September *. IVM. 
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Figure 1 . — (a) Diamond* graphite intergrowth in arti- 
ficially shocked graphite. Note the anhedral mor- 
phology of tbo u tpr diamond grains (D). (b) Typical 
flake from the graphite starting material. The texture 
of this flake may be oompared with that of the un- 
transformed graphite in (a). 


1326 



GENERAL 


1327 



Figure 2. — Unidentified species present in both starting 
maWial and shock product. The cell dimension of these 
cubes is similar to that of the “carbon II” phase re- 
ported by Aust and Drickamer (4.) 

form of euhedral cubes (Fig. 2) which were single 
crystals of uncertain origin. The cell dimension 
of these simple cubic crystals (at, =5.55 A) corres- 
ponds very closely to that of the “carbon II” 
phase (a 0 = 5.545 A) reported by Aust and Drick- 
amer ( 4 ). Inasmuch as the cubes were present 
in very small numbers (< < 1 percent) in both the 
starting material and product we cannot say that 
they were produced by shock or, in fact, that they 
were composed of carbon. It seems very unlikely 
that some extraneous material would possess 
exactly the same cell dimension and crystallo- 
graphic habit as “carbon II”. Yet the cubes’ 
euhedral appearance, the general perfection of 
these crystals, and their low abundance suggest 
that they were contaminants and not shock- 
formed. It could also be that these cubes were, 
indeed, “carbon II” which had grown relatively 
slowly during preparation of the original artificial 
graphite block. Unfortunately, all three alterna- 
tives (contamination, shock-formation, or slow 
growth) seem equally unlikely and it is impossible 
to decide among them with the available data. 

The anhedral morphology of the artificial dia- 
mond might be expected from the short time 
available for crystal growth under the experi- 
mental shock conditions. Shocks generated dur- 
ing asteroidal collision (/) on the other hand, no 
doubt were of 'onger duration and might well be 
expected to result in the formation of subhedral 
dial .ond grains. These grains would probably 
be polycrystanine aggregates, since the shock- 


induced conversion of graphite to diamond ap- 
parently proceeds by a mechanism other than the 
orderly diffusion and addition of carbon atoms to 
diamond nuclei ( 1 , 2 ). v 

X-ray fluorescence analysis of the interior of 
the Canyon Diablo nodule indicated the presence 
of only iron and sulfur. The only carbon phase 
detected was normal graphite. The cubic phase 
seen in the starting material and shock product 
was not observed in the Canyon Diablo nodule. 

The sulfur in the Canyon Diablo material 
apparently exists in two forms; as free rhombic 
sulfur and as troilite (FeS). Figure 3a shows 
anhedral single crystals of sulfur in a graphite 
matrix at an electron accelerating voltage of 
75 kv. At 100 kv, the normal operating voltage 





Figure 3. — (a) Sulfur (S), graphite mixture from Canyon 
Diablo (75-kv electrons), (b) The same area examined 
with 100-kv electrons. The higher accelerating voltage 
has sufficiently increased the temperature to melt the 
sulfur grains. 
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in these studies, the temperature reached was 
high enough to melt the sulfur (m.p. il3°C), and 
further aided in its identification (Fig. 36). The 
troilite in the graphite nodule was identified by 
both x-ray and electron diffraction. It was pre- 
sent in amounts greater thai about 5 percent 
while the sulfur was presen. in lesser amounts 
(about 2 percent). 

The coexistence of troilite and elemental sulfur 
has been observed previously in Type I carbon- 
aceous chondrites ( 5 ) although in these meteorites 
there is considerably more sulfur than troilite. 
Free sulfur has not previously been reported in 
iron meteorites. It may be that the sulfur was 
formed by terrestrial weathering of the troilite, 
the iron thus liberated forming amorphous lim- 
onite (6). An alternative is that these minerals 
are preterrestrial. If this was the case it is 
expected that the S*VS* 4 ratios in the troilit •. id 
sulfur may differ. 
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A.C. CONDUCTIVITY OF A PLASMA* 


C. S. SHENf AND R. L. W. CHEN* 


This paper computes the electrical conductivity of a fully ionized, spatially homogeneous 
plasma under the influence of uniform, periodically alternating electric field. The velocity 
distribution of the electrons is determined by solving the linearized Fokker-Planck equations. 
All the terms in the collision integral are retained, including those representing electron-electron 
interactions. The resultant values of conductivity is expected to be valid in the range of fre- 
quencies from zero to below the plasma frequency. 


I. INTRODUCTION 

The purpose of this paper is to calculate the a.c. 
conductivity of a spatially homogeneous plasma 
using the Fokker-Planck equation. The d.c. 
conductivity of a plasma has been calculated in 
the well-known works of Cohen et al. (1950) and 
Spitzer Jr. and Harm (1953). Their results are 
in good agreement with the later experimental 
works of Lin et al. (1955). Bernstein and 
Trehan compute the a.c. conductivity assuming 
a Lorentz gas model (1960). The a.c. conduc- 
tivity of a real gas should approach that of a Lor- 
entz gas at high frequencies (see detailed discus- 
sions in Section 4). Toward lower frequencies 
their departure is expected to increase so that their 
ratio becomes nearly 2 in the d.c. limit, in accord- 
ance with Cohen et al. and Spitzer and Harm. 
The recent works on a.c. conductivity by Dawson 
and Oberman (1962) consider the time variation 
of the two-particle distribution, which is neces- 
sary when dealing with a.c. currents of ultra-high 
frequencies. However, the domain of applica- 
bility of their work is limited to frequencies much 
higher than the collision frequency. Thus, a 
more precise calculation for the low and inter- 
mediate range of to appears desirable and we pro- 
ceed to do this in accordance with methods to be 
described in the next section. After completion 
of most of the numerical work, a paper by Robin- 
son and Bernstein (1962) came to our attention. 

♦Published in Plasma Physics (Journal of Nuclear Energy Part C), 
6 : 389 - 300 , 1064 . 

tNational Academy o; Science* — National Research Council Re- 
search Associate with NASA and Research Fellow of Institute of 
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They computed the a.c. conductivity using a 
variational technique. Our results obtained by 
direct integration of the Fokker-Planck equation 
will be compared to theirs in Section 4. 

We begin with the Boltzman equation: 


dt r dr m dv ' \ 5t ) c 


d) 


where /, is the distribution function of particles 
of type i, is the change of f { produced by 
collisions. 


Equation (1) is deduced from Liouviiie theo- 
rem to describe a many-particle system under two 
assumptions: 

(i) That the characteristic dimensions of the 
inhomogeneities are much larger than the average 
impact parameter for the particles participating 
in the collision. 

(ii) That the characteristic time variation of the 
process is much longer than the duration of an 
average collision, or in other words, a collision is 
completed and the correlation function is ‘relaxed’ 
before the distribution function itself makes any 
appreciable change. 

It should be noted here that the term ‘duration 
of collision’ is different from the so-called ‘collision 
time’; collision time is the time between two colli- 
sions. For particles interacting through long- 
range forces, this time may be regarded as the time 
in which deflexions gradually deflected the con- 
sidered particle by 90°. Duration of collision is 
the time during which an interaction takes place. 
In a plasma it is of the order to, -1 . In Fig. 1 a 
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> Region one -< 
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0 = 2 Jdv'/,(v')|v-v'| 

(4) 

* COc 


CUp 

M . CD— ► 

> Region two 
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Figure 1. — The ranges of validity computed a.c. conduc- 
tivities. The values of a.c. conductivity obtained in 
this paper is valid in region one. When omega exceeds 
u e — region two — the values calculated by Dawson 
and Oberman begin to be valid. Here o> c is the colli- 
sion frequency, u p is the plasma frequency. 

time-scale diagram is drawn, and the validity of 
our calculation and those of Dawson and Ober- 
man are indicated. 

The explicit expression of 

the nature of the interaction force. In a fully 
ionized plasma, the particles interact through the 
long-range Coulomb forces. The cumulative 
effect of ‘weak’ deflexions resulting from the rela- 
tively distant collisions outweighs the effect of 
occasional large deflexions due to relatively close 
collisions, so one may neglect the contribution by 
those very close encounters (Cohen et al .) — en- 
counters which result in deflexions of 90° or larger. 

Also, the effect of distant particles lying outside 
the Debye length Xp may be neglected because of 
the shielding of inner particles. Thus, in the 

computations of > it is only necessary to 

consider the collisions with impact distance inter- 

* 

mediate between Xp and 6 0 , where b 0 =— — is the 

KT 

impact parameter yielding a 90° deflexion. The 
effects of these collisions are cumulative, and the 
total deflexion produced in an interval of time is 
similar to that of the Brownian motion; hence, 

one may expand i Q powers of [A?], where 

[Av] is the average velocity change due to collisions 
(Cohen et al., Chandrasekhar, 1943). This 
procedure leads to the following Fokker-Planck 
collision integral (Rosenbluth et al., 1957) : 

where 

/ dvy/voiv-vr 1 (3) 


A (A 

V«A 


depends on 


The summation in h t and g sums over all species, 
m t is the mass of the i-th species, e is the electronic 

charge, Xp - ( ) is the Debye length, and 

\tirne 2 / 

v (h is the relative thermal velocity. 

In this paper we consider only plasma with 
singly-charged ions. The extension of the present 
method to those with multiply-charged ions is 
straightforward. 


2. DERIVATION OF EQUATIONS 
AND FORMULAE 

If the distribution function f has an azimuthal 
symmetry about a certain axis, then, following 
Rosenbluth et al., the collision te-m may be 
written down explicitly in spherical polar co- 
ordinates in velocity space : 



where ju*cos 6 is the direction cosine between v 
and E. Equation (6) is an exact expression of the 
Fokker-Planck equation in spherical co-ordinates 
for a distribution function with azimuthal 
symmetry. 
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Wc assume that the system is subject to a weak 
electric field Eoe*" 1 whose direction lies along 
z-axis. Then following Chapman and Cowling 
(1939) and Spitzer Jr. we expand /, in a power 
series of E: 

Mr, t ) -/,«(*) +1 W\y, t) *) + •••( 7) 


^(x)=/e~ , ' , di/ 

o 

(15) 

tf(x)=*(x)-x*'(x) 

(16) 

X 

7»(x) = / y n D(y)e~ v 'dy 

(17) 


0 


where / 1 <0) (») is a time-independent Maxwellian 

distribution and f (I, (r, t), / e> (v, t), are the 

perturbed part due to applied electric field. When 
a steady state has been reached and no transient 
current exists, the time-dependent part of f, U) 
(v, () must be proportional to e*" 1 . Since the 
average energy imparted to the electrons between 
encounters is small compared with their kinetic 
energy, the velocity-dependent part of /< <fl (v, t) 
can be written as e~ lm * PI2k,1 D t ll) (v)y. Therefore, 
we have 

f( Av, t ) = ( 2 ^Sji 7 ( 8 ) 


Combining equations (1), (6) and (7) keeping 
only terms linear in Et>, we obtain 


a (1 > - 

Q/< j C fa 

dt 'm x 



where ( 6f i /6i) c , (1> is the linearized Fokker-PIanck 
collision integral. 

Since the ions’ contribution to electric current 
is negligible compared to electrons, we will con- 
sider only electron distributions and drop the 
subscript i in the distribution function hereafter. 

Substituting equation (8) into equation (9) 
we find, after some algebraic manipulations, the 
following second-order linear integral-differential 
equation: 


D"(x) +P(x)D'( x ) +Q(x)D( x ) = «(x) +S( x) (10) 
where 


p / v ) _ _ gy — - | — (*) 

{x) X x + H(x) 


Q(x) = 
R(x) = 


— iB\ +2(l-f$=2x $*) . 1 

H(x) x ! 

2d x 4 8(2.4x‘-2 x 4 ) 


H( x ) 3t 1/1 H (x) 


16 


( 11 ) 

( 12 ) 


/«(») (1«) 


S(x) t x/»(x) - 1 *2 xh(x) 

-X 4 /,(x)(1-1-2x*)} (14) 


Io ( 00 ) = jD(y)e v 'dy 
0 

with 

a = — EoK T/iren In X X = Xj>/6o 



4'%/2Xw 
u p In X 


(18) 


u 



where u e is approximately the 90° deflexion time 
of a particle with thermal velocity. When u = 0, 
equation (10)* reduces to equation (8) of Spitzer 
Jr. which considers d.c. electric conductivity. 

In a d.c. electric field, the electrons are not 
accelerated in a steady state. Hence, the inertia 


q } 

force term is zero and 7 0 ( « ) = _ a . 

8 

3. SOLUTION OF EQUATION 


Equation (10) is a linear integral-differential 
equation whose unknown D(x) is a complex func- 
tion of a real variable. The present section will 
discuss the method of its solution. As will be 
evident in what follows, the procedure for 
numerical integration is far from straightforward. 

On the one hand, we encounter the problem of 
the instability of the solution at small and at 
large x- Because of the existence of singularities 
in equation (10) at x = 0 and at x= °°, a slight 
deviation of D(x) at either small or large x, tends 
to be built up quite rapidly. In order to obtain 
a physically acceptable solution, it is required 
that D(\) does not a~oroach infinity too fast, lead- 
ing to infinite conductivities. The starting value 
of D at small x can be obtained by means of a 
series solution. Because of the instability, we 
cannot proceed to integrate in a step-wise manner. 
To overcome this difficulty, we adopted a scheme 


♦Not© that /♦(•) it es s ential ly the total change of momentum of 
electrons arising from eiectronion interactions. Since the mutual 
electronic interaction cannot change the total momentum of the 
electrons» /•(•), by Newton's second law, must equal the total force 
exerted on the electrons by the applied held minus the inertia force of 
electrons. This relation gives us: 

* , . 3ir* t B r . . 

g /*(*). 
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used by Cohen et al. We shall refer to their 
paper for full details. 

On the other hand, we note that 7 0 (co) is no 
longer a known quantity as it is in the case ox d.c. 
conductivity; it depends on the solution D(x) 
itself. We proceed as follows: Since we w;i.nt 
conductivities at different frequencies, it is neces- 
sary to obtain solutions for different values of the 
parameter 5. We begin with a small value 
5 = 0*05. Using an 7 0 (<») taken from the d.c. 
case, i.e. / 0 («) = 0*655, we obtain a solution to 
equation ( 10 ) from which we get a new 7o(«). 
Next, we pass on to 5=0*1 using the io(°°) ob- 
tained for the previous 5. In this way, we pro- 
ceed to ever-increasing values of 5, until the initial 
adopted 7 0 (® ) and the final calculated 7 0 (o°) 
differ by no more than 2 percent. This occurs at 
5 = 1*37. From this point on, we resort to a 
method xf systematic trials. The initial and 
final /„(«>) for all values of 5 agree to within 2 
percent, which is considered sufficiently accurate 
for the present purposes. 


4. RESULTS AND DISCUSSION 

The current is given by 

J= — eJ"dw/j[(v, t) = AEoe <a, Ii(co) (19) 


where 


/»(»)= Jx*7>(x)e X dx 
, _2 (2KT) m 

A —ft 


3 * a '*m l '*e s ln X’ 

Since J =<rE, we have the complex conductivity 
ff<=AI t (c o), 

the impedance 

Z __J £_ 

A7,(co) /,(»)’ 

the resistance 


( 20 ) 

ivity 

( 21 ) 

( 22 ) 


r, 1 Relt(a>) cReh(a>) 
A |/ J ("o)|* = |I,(»)|” 


(23) 


and the reactance 


with 


1 7m7,(<x.)_c7m/,(®) 

X A TKC^TT 5 IT^yp (24) 


1 3 1 

c — r=« 


A 2 (2KT) iia 


We* In X. 


\ 


It may be remarked here that the a.c. conduc- ^ 
tivity depends on three factors: 

(i) The inertia of the conducting electrdns. 

(ii) The mutual interaction among electrons 
and ions. 

(iii) The mutual interaction among electrons 
themselves. 

The mutual electronic interactions have no 
direct effect orx conductivity since the total change 
of momertum due to such interactions is zero. 
Nevertheless, they alter the distribution of elec- 
trons and thereby modify the effect which elec- 
tron-ion collisions and electron inertia have in 
impeding the current. When a is small, the con- 
ductivity is primarily determined by collisions. 
The inclusion of electron-electron interactions 
reduces the conductivity by a factor of approxi- 
mately two. As a increases and becomes of 
order u e , this effect becomes less and less impor- 
tant because there is then insufficient time in each 
a.c. cycle to allow an effective modification of the 
. jtribution by electron-electron interactions. 
When u well exceeds w e , we may neglect this 
effect and D(x) reduces to 


7><(x) = 


. , iB 3 
i+yx 


(25) 


and the corresponding conductivity becomes 

7 7 -* 1 

a,ix)=A j JL ih ix (26) 


» l+yx- 


which is just the a.c. conductivity of a Lorentz 
gas (Bernstein and Trehan). 
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If we further increase «, the inertia of electrons 
become dominant. Then we may treat collision 
effect as a perturbation and obtain 



(27) 

(28) 

(29) 

(30) 


In Table 1 the values of D(x) for B= 1 are 
given and compared with the corresponding 
values for the d.c. case obtained by Spitzer Jr. 
and Harm. 

In Table 2 the resistance, the reactance and the 
absolute value of conductivity are given for vari- 
ous B from 0 to 10. For 2? >10, one may use 
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Table 1 . — Valuesof D(x)Joro>—v t andu>— 


X 

03 


W ““0 


Re ID (*)1 

Im[D{x) ) 

D(x) 

0. 10 

0.0005887 

- 0.0002029 

0.0008093 

0, 11 

0. 0009252 

- 0.0003438 

0.001300 

0. 12 

0.001376 

- 0.000542 

0.001970 

0. 13 

0.001956 

- 0.0008068 

0.002847 

0. 14 

0.00268 

- 0.001149 

0.003955 

0. 15 

0.00356 

- 0.00158 

0.005317 

0. 16 

0.00461 

- 0.00210 

0.006955 

0.17 

0.00583 

- 0.00272 

0.008886 

0. 18 

0.00724 

- 0.00346 

0.01113 

0. 19 

0.00884 

- 0.00431 

0.01370 

0.20 

0.01063 

- 0.00528 

0.01660 

0. 22 

0.01483 

- 0.00761 

0.02347 

0.24 

0.01985 

- 0.01048 

0.03180 

0 , 26 .... 

0.0257 

- 0.0139 

0.04165 

0.28 

0.0324 

- 0.0180 

0.05304 

0.30 

0.0400 

- 0.0226 

0.06601 

0.32 

0.0483 

- 0.0279 

0.08057 

0. 34 j 

0.0575 

- 0.0339 

0.09672 

0.36 ! 

0.0675 

- 0.0405 

0.1145 

0. 38 

0.0783 

- 0.0478 

0. 1338 

0. 40 

0.0899 

- 0.0557 

0. 1548 

0.44 1 

0.1153 

- 0.07366 

0.2015 

0.48 

0. 1435 

- 0. 09436 

0.2545 

0.52 i 

0. 1744 

- 0. 1179 

0.3137 

0. 56 

0.2080 

- 0. 1442 

0. 3792 

0.60 

0.2439 

- 0. 1734 

0.4508 

0.64 

0.2882 

- 0.2055 

0.5285 

0.68 

0. 3227 

- 0. 2405 

0.6123 

0.72 ... 

0.3652 

- 0. 2785 

0.7023 

0.76 

0. 4096 

- 0.3196 

0.7983 

0.80 | 

0.4559 

- 0. 3637 

0.9005 

0.88 

0. 5535 

- 0.4813 

1.123 

0. 96 

0.6570 

- 0. 5718 

1.371 

1.04 

0.7656 

- 0. 6957 

1.645 

1.12 

0. 8782 1 

- 0. 8335 

1.945 

1.20 

0.9937 

- 0.9858 

2.273 

1.28 ! 

1.111 

- 1.1531 

2.630 

1.36 

1.2290 

- 1.3359 

3.017 

1.44 

1.3457 

- 1.5347 

3.435 

1.52 

1.4598 

- 1.7500 

3. 887 

1.60 

1.5693 

- 1.9820 

4.375 

1.76 

1.7657 

- 2.4962 

5.465 

1 * 92 . ... 

1.1915 

- 3.0739 

6.728 

2.08 

1.9973 

- 3.7049 

8. 190 

2. 24 

1.9913 

- 4.3694 

9.880 

2.40 

1.8852 

- 5.0382 

11.83 

2.72 

1.4266 

— 6 . 2*97 

16.62 

2. 88 

1.2198 

- 6.6809 

19.53 

3.04 

1.3105 

- 6.8066 

22.74 

3. 20 

2.1113 

— 5.6758 

26.00 
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Table 2. — The conductivity, the resistance and the reactance of a.c. current 



<r/A 

R/c 

X/e 


Real gas 

Lorentz gas 

Real gas 

Lornetz gas 

Real gas 

Lorentz gas 

0,0... 

1.734 

3.0 

0.577 

0. 333 

0.0 

0.0 

0.05-. 

1.729 

2.880 

0.577 

0.340 

0.045 

0,70 

0.1— 

1.713 

2.653 

0.577 

0.354 

0.089 

0. 130 

0.15.. 

1.687 

2.430 

0, 578 

0.367 

0. 134 

0. 185 

0.2... 

1.651 

2.233 

0.579 

0.380 

0. 178 

0.236 

0.25.. 

1.608 

2.061 

0.58 r 

0.393 

0.223 

0.285 

0.3... 

1.561 

1.913 

0.582 

0.404 

0. 267 

0. 332 

0. 35. _ 

1.510 

1.784 

0. 584 

0.425 

0.311 

0.377 

0.4... 

1.458 

1.671 

0.587 

0.425 

0.354 

0.421 

0,45.. 

1.406 

1.572 

0.590 

0.434 

0.397 

0.465 

0.5... 

1.354 

1.483 

0.593 

0.443 

0.440 

0.508 

0.55.. 

1.303 

1.405 

0.597 

0.452 

0. 482 

0.550 

0.6... 

1.255 

1.334 

0.600 

0.460 

0. 525 

0.592 

0.65.. 

1.208 

1.270 

0.603 

0.467 

0.567 

0, 634 

0.7... 

1.164 

1.212 

0.607 

0.475 

0.608 

0.675 

0. 75. . 

1.122 

1. 159 

0.610 

0.482 

0.650 

0.716 

0.8... 

1.082 

1.111 

0.614 

0.488 

0.691 

0.756 

0.85.. 

1.045 

1.066 

0.617 

0.495 

0.732 

0.797 

0.9. . 

1.009 

1.025 

0.620 

0.501 

0.773 

0.807 

0. 95.. 

0.976 

0.987 

0.624 

0.507 

0.813 

0.877 

1.0— 

0.944 

0.952 

0.627 

0.513 

0.853 

0.916 

1.1— 

0.886 

0.889 

0.634 

0.524 

0.934 

0.996 

1.2 

0.834 

0.835 

0.640 

0.534 

1.013 

1.074 

1.3 — 

0.786 

0.785 

0.645 

0.544 

1.094 

1.159 

1.4 

0.743 

0.741 

0.750 

0.553 

1.172 

1.230 

1.5— 

0.721 

0.702 

0.657 

0.561 

1.222 

1.308 

2.0— 

0.560 

0.557 

0.672 

0.599 

1.588 

1.693 

3.0--. 

0.408 

0.398 

0.711 

0.655 

2.335 

2.453 

4.0-.. 

0.313 

0.305 

0.734 

0. 696 

3.086 

3.207 

5.0— 

0.255 

0.248 

0.766 

0.728 

3.842 

3.960 

6.0— 

0.214 

0.211 

0.784 

0.754 

4.601 

4.710 

7.0 — 

0. 183 

0. 182 

0.790 

0.776 

5.380 

5.460 

8. 0 

0. 162 

0. 161 

0.796 

0.794 

6.08 

6.21 

9.0 — 

0.144 

0.143 

0.814 

0.810 

6.69 

6.81 

10.0... 

0. 130 

0. 129 

0.830 

0.824 

7.51 

7.6 


equation (26) to compufu them. The error will 
be within 2 percent. For B>50' the collisions 
become unimportant and equations (27)— (30) will 
give the correct values to within 2 percent. How- 
ever, there the validity of the Fokker-Planck 
equation already becomes questionable and one 
should use Dawson-Oberman’s values instead 
of ours. 

In Table 3, the complex conductivity calculated 
in this paper are compared with those obtained 
by Robinson and Bernstein. They showed 


that transport coefficients obtained from the 
Fokker-Planck equation should possess an ex- 
tremal nature, and proceed to calculate conduc- 
tivities using the variational technique. T heir 
Table 6 gives conductivities for various values of 
the logarithm of &>/<•>< including very large values 
of u e . In our Table 3 only those values are in- 
cluded for comparison which fall within the range 
of validity of the Fokker-Planck equation. The 
discrepancy between the two results are generally 
within 5 percent. 
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Table 3. — Comparison of the a.c. conductivities obtained in this paper ( direct integration) 
with those obtained by Berstein and Robinson ( variational calculation) 


w/wc 

Log X 

(X-4V2W/W) 

Rea/ A 

—lma/A 

Bernstein 

Robinson 

Shen-Chen 

Bernstein- 

Robinson 

Shen-Chen 

0.0. 

— CO 

1.734 

1.734 

0.000 

0.000 

0.0057 

-3.0 

1.734 

1.734 

0.004 

0.002 

0.0179 

-2.5 

1.730 

1.732 

0.G49 

0.408 

0.0565 

-2.0 

1.705 

1.722 

0. 154 

0. 151 

0. 179 

-1.5 

1.605 

1.608 

0.447 

0.443 

0.565 

-1.0 

0.992 

0.987 

0.816 

0. 813 

0,901 

-0.8 

0.637 

0.632 

0.785 

0.787 

1.425 

-0.6 

0.357 

0.353 

0.545 

0.649 

2.261 

-0.4 

0. 183 

0.179 

0.477 

0.482 

3.559 

-0.2 

0.087 

0.084 

0.329 

0.336 

5.65 

0.0 

j 0.040 

0.038 

0.201 

0.204 


We should add that our results can be readily 
applied to the case of conductivity in the presence 
of a uniform magnetic field. The addition of the 
magnetic field leads to equations which are en- 
tirely similar to (10). If the electric field is 
parallel to the magnetic field, the conductivity is 
not affected. If it is perpendicular to the mag- 
netic field, the conductivity becomes a H {o>) m <r 
(a+'w), where <r is the function obtained in this 
paper and ug^eH /me. 
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PARTIAL WAVE THEORY OF DIATOMIC MOLECULES, 

II* 

, A. TEMKIN AND A. K. BHATIAt 

The method of expanding the wave function of diatomic molecules in a series of orbital angular 
momentum eigenfunctions (partial waves) h extended to two-electron homonuclear molecules. 

The angular momentum is a function of the Euler angles only; a symmetric cl Ace of those angles 
is used which greatly facilitates the description of the exchange character of the wave function. 

As a result, explicit equations for the 3 dimensional “radial” functions can be derived for all the 
different magnetic parity, and exchange states. For £ states t.io partial wave sums only go 
over alternate values of the angular momentum giving rise to the expectation of an even more 
rapidly convergent series in these cases. 


I. INTRODUCTION 

In a previous paper 1 as well as others,* the idea 
of expanding the wave function of a (homonu- 
clear) diatomic molecule in a series of angular 
momentum eigenfunctions was introduced for the 
purpose of affording a very natural method of 
successive approximations in these problems. 

Previous calculations 1 * using this idea have 
almost exclusively been confined to the one-elec- 
tron molecule H* + . In that case, the formalism 
is rather easily constructed, and the results are in 
the nature of a check as the problem is amenable 
to exact calculation . 1 The comparison with the 
exact results is nevertheless most encours'ing . 1 

The main immediate goal of this approach is 
two-electron molecules for which there is no 
analytically exact solution. As a prerequisite for 
effecting this approach, however, it is necessary 
that two electron atoms be completely understood. 
Although the general approach to the two- 
electron atomic problem is quite old 4 , only results 
for specific angular momentum states have until 
recently en derived. The difficulty lay in the 
description of the character of the wave function 
under the exchange of electron coordinates which 
becomes very complicated 4 with the original 
choice of Euler angles*. The explicitly complete 
solution has required a symmetrical choice of 


‘Published as Goddard Spec* Flight Ctnior Docwnt X *6404+46, 
March 1MH. 

tGoddard Space Flight Center; National Science Foundation, Na- 
tional Research Council resident research associate. 


Euler angles, which Holmberg* has provided, and 
has therefore only just been carried out 7 . 

In section II, we very briefly review the main 
results of that investigation and present some 
additional transformation properties of the Euler 
angles and vector spherical harmonics which bee* 
on the present application. In section III we 
describe the construction of the molecular wave 
functions, and in section IV we derive the radial 
equations. These equations couple states of dif- 
ferent angular momentum through the potential 
energy, but the kinetic energy acts in the same 
way as in the atomic case. The new feature is, 
therefore, only the coupling terms of the molecu- 
lar potential and this in turn depends on integrals 
over three vector spherical harmonics. The heart 
of section III is the resulting formula for these 
integrals and the selection rules for its non-vanish- 
ing. This latter consideration leads to the 
decoupling of the Z states accordirg to even (+) 
or odd (-)l. This decoupling is not present for 
other m states corresponding to the absence of the 
± quantum numbers in those cases; thus the Z 
states provide particularly favorable cuses for the 
rapid convergence of the partial wave expansion. 

Section V presents the zeroth-order equation 
for the states and discusses a specific calcula- 
tion which can be related to the accuracy of this 
approximation and which by inference bears upon 
the convcrgenr properties the entire sequence 
of approximations which thi* paper projects. 
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II. THE EULER ANGLES AND 
TRANSFORMATION PROPERTIES 

The Schrodinger equation of a two-electron 
homonuclear molecule whose nuclei are con- 
sidered fixed (Bom-Oppenheimer approximation) 
will deper.d on the two vectors of the electrons Ti 
and I* which can be measured relative to an origin 
defined an the midpoint between the nuclei and 
parametrically on the intemuclear separation R.<o. 

In place of the six coordinates ri, rj one can in- 
troduce the three Euler angles and three residual 
variables; the Euler angles describe the orienta- 
tion of the plane containing the two unit vectors fi 
an*i f» relative to a space-fixed coordinate system 
(whose z-axis is by definition the intemuclear axis) . 

The particular Eu’er angles we use where intre- 
% duced by Holmberg® and have the advantage of 
being symmetrically disposed with respect to the 



Fioohb I. — Perspective drawing of (Holmberg’s) Euler 
angles and the unit vectors of the problem. 


tw>electron. Figure 1 contains a diagrammatic 
'^presentation of these angles and the unit vectors 
with which we shall be concerned. 

The axes are defired by 


«“(fiXfj)/sin 9u 

(2.1) 

H ■ (iXiO/sin 9 

(2.2) 

f-VXi! 

(2.3) 


where 


The Euler angles are obviously: 0 = angle between 
* and 4* = angle between i and x'; ^ = angle be- 

tween i' and (ri— ?*)- Trigonometric relations 
between these angles and the particles spherical 
angles have been worked out in reference 7; the 
reader is encouraged to familiark * himself with 
that discussion, although we shall attempt to 
make this treatment reasonably self-contained. 
Specifically, we shall repeat he r < 1 the relation of 
the x, y and z components of fj and r* to the Euler 
angles in order that the transform m properties 
given below can readily be checked. 

— =sin cos ^i = cos <t> sin (^— x A6\i) 

T 1 

+cos 9 sin -f> cos (v/ — x A9\t) (2.5) 

— = sin 0i sin #>i = sin $ sin (4 / —' l A6u) 

—cos 9 cos 4> cos (^— H^u) (2.6) 

— = cos 0i® — sin e cos (p — J40u (2.7) 

T i 

The formulae for analogous components of fj may 
be obtained from the above by letting 
The angle is the angle between f t and ?*. It is 
not one of the Euler angles, but rather one choice 
of residual coordinate, tl : other two usually being 
chosen as the two distances ri and r s . The resid- 
ua! coordinates are often called generalized radial 
coordinates; an alternate choice is ri, „ r^, the 
letter being the inte r ele<vron disiinc.e. 

The distinguishing feature of the Luler angles 
as opposed to residual coordinates is the fact that 
the angular momentum depends only on the Euler 
angles. In this connection, it is well known that 
square of the total orbital angular momentum, 
M’, can be conserved in atoms but not in diatomic 
molecules, whereas the z-componont of the angular 
momentum, M„ is a possible constant of the 
motion in both cases. A cent: ..1 role in the atomic 
problem is therefore played by the eigenfunctions 
of the total angular mor.vnvum Dj M (0,<fr,'k): 

*, i) ( o ) ( 2 . 8 ) 
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0= \k+m\ + \k-m\ +2 (2.10) 

and 

w m 1 r (2/-M) (*H|H-m|+lH-m|)t v (/-l|H-m|H-ij*-ml)n » 

,nk ( | Ac — m|) !|_ 8 t* (^-i|A+m|-§|*-m|)! X (/+§|*+m|-i|ife-mi)!j ^ ll) 


The eigenfunction character of these functions 
is then given by 

M hH(t + 1 )D-* (2.12) 

where the explicit form of M l and M, has been 
derived in reference 7 and in any case is known* 
In the atomic problem if we are describing a state 
of definite orbital angular momentum l with z- 
component to, then the wave function is at most a 
finite sum over the (2£+l) function for the given 
t and to: 

'r,„(r,,r,) = 2 (2.14) 

»*-< 

Actually the above sum is more general than 
what is needed. Specifically one has an addi- 
tional operation which commutes with the Ham- 
iltonian and therefore gives rise to an additional 
constant of the motion. This is the operation of 
space inversion (called parity, p, in the ator ..o 
application 7 ), usually labelled iE in mole alar 
applications* corresponding to the operation fi 
— Ti r s — ►— r* in all cases. It turns out that 
only tue Euler angles are involved in this trans- 
formation, and if we denote them by a three 
component vector, (0, 4>, 4'), then 7 

iE(*\ = ( 4> \ (2.15) 

W W*/ 

and 

« (- l)*Z) < w *(ff,$,'F) (2.16) 

Therefore, the wave function (2.14), if it is to be 
an eigenfunction of iE, must be restricted to even 
or odd k . Since the action of the Schrodinger 
equation is to couple all possible gf-* functions 
for the given l and to, it can be seen that the in- 
vocation of parity halves the number of functions, 
and thus among other things is of considerable 
practical utility. 


In the case of diatomic molecules the invariance 
against rotations, except those around the i-axis, 
is lost. Thus according to the above arguments 
only the rotation quantum number to survives. 
The Idea ef the partial wave theory, then, iete 
expand the molecular wave function in terms of 
the orbital angular momentum: 

t-UI 

The above sum too may be delimited according 
to the symmetry properties of the states we intend 
to describe. These symmetry operations are well 
known for the diatomic molecule,*, and except 
for exchange, ««, they all involve only the Euler 
angles. The transformation involved in the 
residual variables for tn, is r^r* independent of 
whether the third variable is taken as On or r«. 

In Table I we have summarized the transforma- 
tion properties of the Euler angles and the vector 
spherical harmonics under some of the major 
symmetry operations of the two-electron homo- 
nuclear diatomic molecule. 

In the next section we shall utilize these proper- 
ties to construct wave functions with the desired 
symmetry properties. We emphasize that not 
all these properties are independent. Thus for 
example a full space inversion is the product of 
the inversion about each of the axes separately: 

iE = <r*<r px& ty 

where the factors on the right can be written in 
any order. Therefore, the operation <r h is 
redundant. 

The effect of additional symmetry operations 
such as <r, on the angular momentum eigenfunc- 
tions and variables can easily be worked out from 
Table I by noting, in this case, that 

io,=iE’<r, 

Not all symmetry operations commute with M,. 
Specifically <r, does not commute with M,. This 
may be easily ascertained by envisaging a vector 
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Table 1. — Symmetry properties of the Euler angles and the angular momentum eigenfunctions. 


Transformation Property 

e 

* 

* 

V 

iE space inversion (parity) 

e 

♦ 

1 

*+x 

(- 

<Tk reflection through x = y plane 
(inversion of z-coordiaates 

e 

♦ +X 

*+ir 

(- i)* + “z>7-* 

Or* reflection through z — x plane 
(inversion of y-coordinates) 

tc-q 

2x— ♦ 

♦ 

( — 

Ow . reflection through z — y plane 
(inversion of x-coordinates) 

T-0 

X— ♦ 

♦ 

(-1 ) <+ — *i>7“-* 

i 

cii exchange of particles* 

T-0 

♦ 4-x 

2 T-* 



*Thi« operation also affects the residual coordinates. See text. 


whose projection in the x-y plane is in the first 
quadrant. If first you reflect it through the y-z 
plane thus putting in the second quadrant 
and then rotate around the 2 -axis by 90°, this 
vector will end up with a projection in the third 
quadrant. If first you rotate by 90°, putting it in 
the second quadrant, and then reflect through the 
y-z plane, you put it back into the first quadrant. 
Thus 

[<r„ M.1^0 (2.19a) 

This means that in general one can’t construct 
simultaneous eigenfunctions of M, and ay 
There is one exception to that situation, and 
that is 2 states. 

For in that special case 

M,D°*=0, (2.1V) 

so that 

0, (2.18) 

or 

k„ itf.] w -o=0 (2.19b) 

We shall see in the next section that because of 
this fact the structure of the 2 wave function and 
the resulting 2 radial equations is quite different 
from all other states. 

III. WAVE FUNCTIONS AND THE POTENTIAL 
ENERGY 

In order to facilitate the derivation of the radial 
equations, it is convenient to construct an alter- 
nate set of orthonormal vector spherical har- 


monics which for m = 0 are real 7 . We define 

D}" , ' ,,+ =(V2+j- 0 ,( 2-V2))" 1 

(D?"+D?-) (3.1a) 

D l r- x) ~=(V2ir i (D7- K -D?--‘) (3.1b) 

The important point about these alternate sets 
of vector spherical harmonics is that they are 
eigenfunctions of exchange: 

tl2 D^-«>= t =±(-l)‘D<»-*)± (3.2) 

while still retaining their eigenfunction character 
with respect to parity. M 2 , and M, (with eigen- 
values (—1)*, ft 2 /(£+l), h m respectively). 

For we have already indicated in the 

previous section that the correct expansion of the 
complete wave function is 

S(/i m, * ) »z>r , * >+ 

t-\m\ x 

+ y* m .« , -( r ) D <m,. ) -) (33a) 

The double prime on the summation indicates that 
every second value of the summation index (in 
this case x) is to be taken. This insures as stated, 
that the wave function is an eigenfunction of iE. 
It is precisely this property which distinguishes 
the homonuclear from the heteronuclear mole- 
cule. Thus there is established in the present 
homonuclear case two classes of solutions for a 
given m, one for even and one for odd x. These 
are of course the gerade and unger&de states 
respectively. 
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The radial functions are written as 

functions of r which is short for the generalized 
radial coordinates, either r h i 2l O 12 or r : , r 2 , r n . 
The superscript indicates that the radial functions 
depend on the magnetic quantum number m. 
How" er the well-known degeneracy of the energy 
with respect to ± m indicates that the radial func- 
tions for ±m states must be simply related to each 
other. The relation is undefined up to a ± sign, 
but a consistent choice is (see below) 


fr m - x)± {n,T2M = (-l)7< ( "- ,,± (r 1 ,r l ,M 


The derivation of radial equations divides itself 
into two parts, the kinetic energy and the poten- 
tial energy. The kinetic energy is diagonal with 
respect to l and m but by no means trivial with 
respect to how mixes the allowable x components. 
However, this aspect of its behaviour has beer, 
derived in detail in reference 7 and will not change 
in this application. 

The potential energy in the atomic problem 
depends only on the generalized radial variables 
Ti, r 2 , r u and therefore acts as additional, com- 
pletely diagonal terms in the atomic radial equa- 
tions. On the other hand the molecular potential 
energy is precisely what destroys the spherical 
symmetry of the Hamiltonian. Specifically for 
nuclei of charge z (energy in units of rydbeigs) 
to complete molecular potential is 


.2 2z 2s 2g 2z . 2 

r 12 12.11 Ra2 R«1 Rb 2 Ti2 


(3.4) 


The last term is the interelectron repulsion and 
since it commutes with the angular momentum as 
discussed above need not be considered further. 
For the rest we can expand 


F m w= — 4«2 Q\{?Rab, ri)Px(di) 

X even 

— 4«S OkHRab, rj)Px(ds) 

X even 

where 

i*X 


g*(x, y) = 


x« *<y 


Ii a+i 


y<x 


(3.5) 


(3.6) 


Rab is the intemuclear separation. The Legen- 
dre polynomials P\(d) (note what we mean is 
eg. Pi(d) = cosd not Pi(d)=d) can further be 


expanded in terms of the Euler angles and On 7 : 


Px( * l)= (2xTTpS K + (*i *)D?' M)+ 

It even 

-«r(0i2)l>i o, '‘ ) "] (3.7a) 

Applying « w , we get from the above 

If 

Px(t, 2 )= ( 2 xTTpS 

it even 

-<*r(ffi 2 )Z>f M,_ ] (3.7b) 

where 

F(I-Som) cos(| 0| (3.8a) 

oT(«-(-l)'K(./2)[|^|] } 8m (|0„) 

(3.8b) 

The associated Legendre polynomials are under- 
stood to have the phase of Magnus and Ober- 
hettinger 10 , and their argument is understood to 
be inserted in the transcendental form of these 
functions as parenthesized above. 

Since P\(<f) has parity (— l) x and from (3.5) X 
is even, it follows the sums in (3.7) go over only 
even ft as indicated. 

The point of these expansions is that when the 
wave function (3.3) is substituted into the 
Schrodinger equation 

H* m = E* m (3.9) 

where 

H=T+ V m <st + 2/r « (3.10) 

and T is the kinetic energy, only V m0 < will couple 
terms of different l. The decompositions, Eqs. 
(3.7), then tell us that the coupling from a state 
(radial function) with quantum members t, x, m 
to one with quantum numbers t f , x r , m (the molec- 
ular interaction is still diagonal in m) will be gov- 
erned by integrals over three vector spherical 
harmonics. Two of these angular functions come 
from the coupled states and one from the interac- 
tion; the integral in question is therefore a linear 
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combination of integrals of the form 

n 2t 2r 

J f j (DT/ xi )*Dl-W , t sin OdOdm- ~ (f < X (AX - x* 1 x r ) (3.11) 

0 0 0 f 


The formula for this integral given on the rhs 
of the above equation is well known 11 in terms of 
the Clebsch-Gordan (C-G) coefficients. The se- 
lections for the vanishing of this integral can 
therefore be deduced from those of the C-G 
coefficients. 

Consider first 2 states (m= 0). In this case 
the first C-G coefficient becomes ( t X 00/£'0). 
This coefficient vanishes unless /+X4-£'=even 
(integer). But X is an even integer, Eq. (3.5, 
therefore: 

(+(' - even 

This says that for 2 states only radial functions 
of even l will be coupled together and only odd l 
states will be coupled but there will be no inter- 
mixing. Thus for 2 states the expansion of the 
wave function, Eq. (3.3a), separates into two 
classes defined by the evenness or oddness of l\ 

*o= 22(/} 0, * )+ (r)Z)) 0 -* )+ 

+/< M) '(r)Dj M> ~) (3.12) 

The molecular quantum numbers ± are deter- 
mined by the parity of ( which is the eigenvalue 
of i <j, y . This then constitutes the exception to 
the noncommutivity of <r f „ r. 'tit M, which was 
noted in Eq. (2.19). 

The quantum lables of the 2 states are com- 
piled in the entries of the following simple array: 



For raj^O, the relevant C-G coefficient is not 
necessarily zero when /+X-K' = odd. Thus the 
decoupling according to even or odd l (± sym- 
metry) does not occur, corresponding to the non- 
commutivity of o-, with M„ Eq. (2.19b). Here, 
however, there exists a well known degeneracy in 


energy with respect to ±m eigenfunctions. If 
the eigenfunction for +m is given by (3.3a) and 
that for —m is given by 

tf 

4'_ m = 2 s(/ < ( - m) *> + I> < ( -’ B '*> + 

<-|m| x 

+/ < <- m .‘>-ZV-’".*)-), (3.3b) 

the obvious question is what is the relation of the 
/ / (m .* )± to fc l ~ m ' x)± . The answer to this question 
still resides in the same C-G coefficient, in this 
case in the symmetry property 11 : 

(CkmO \ Cm) = ( — 1 ) <+x+< (Ck--m0\t'—m) (3.13) 

Since X = ?ven, we see the si 0 Ti of the coupling of 
an l to ar. V state is reversed for ^_ m relative to 
what it is for if £-K' = odd (i.e., I and V have 
opposite parity). If the radial functions change 
sign in the same way, then the equations (and 
hence the eigenvalues) will in fact not change, 
i.e., we must have 

(_i)^(m,«)± ( 3 .l4a) 

This relation could be altered by a minus sign; 
the relation 


(_ i)< + »^(»,*)± (3.14b) 

is equally acceptable and has the advantage that 
the first radial function that appears in the expan- 
sion of will be the one which doesn’t change 
sign when used in 

In all cases we have by virtue of the very sim- 
ple exchange property of ZV m '* )± , Eq. (3.2), that 
the complete spatial wave function will be 

symmetric (upper sign) or antisymmetric (lower 
sign) if the radial functions have the exchange 
property 

ft;}) 

= ±(- l)VT** ,+ (ri, ( 3 - 15a ) 







1342 

and 

/r*’-(r„ r, ,{»;;}) 


PUBLICATIONS OF GSFC, 1964: I. SPACE SCIENCES 

IV. RADIAL EQUATIONS 


-±(— 1)' + W— »-(n, r,, {*■’}) 


The foregoing discussion together with the 
(3.15b) knowledge of the kinetic energy operator 7 essen- 
tially allows us to write down the “radial” equa- 
tions for arbitrary m. Specifically substituting 
These properties essentially halve the inde- (3.12) or (3.3a) into the Schrodinger equation, 
pendent variable space over which the equations multiplying on the left by and integrat- 

need to be integrated 7 . ing over the Euler angles, one obtains 

'/(/+!)— x 2 , x 2 eo+fa , 


2 sin 2 6)2 ' 4 2 sin fa 


-2^)j/j w '*> + , 


+ (x+l)(x+2 )^ e *-B t x+ . 2 ft m,x+ 1 )+ (r)+ ~ (l±3go.«-2) A\_ 2 ALn m , x _ 2)+ l 

sm fa W+ 4 sinfa (l + 3 0 .,- 2 (V2-l)) Vu U (r) J 

A cos fa a 1 

\ r * ryi \2 sin (Jfa) 12 2 sin fa/* ^ * <?fa 


5l,(l— 2Al) Am,*)- , (x+l)(x + 2) /(*,«+!)-/ \ 

2 sin fa f ‘ sin fa W 


1 


(1~ V*-2) 


4 sin fa (1+5 o,,-2(v^-1)) 
2 m 


(r) J 


2 2 9 / v 

2 2 9*\o R ab> r /)(l — 252^6, x _)ax* A (fa)/m(^x-(-, y x , 7/)/f W ' v>v (r)=0 

1 f J-l 7 X(even) V 6 / 


(4.1a) 


l(M-l) -X 2 . X 2 cot fa / \ 

sin 2 fa + 4 + 2 sin fa 5l * (1 lA *)y< W 






»>+M ■ SjAV-2Ai) f(mtK)+ 


dd 12 


2 sin fa 


-/r* )+ (r) 


— ^ ^■«+t)(* , + l)(x+2) R Xm, «+2)+/ A, 1 (1 +3i 0 , < — 2 ) ,(m.*-2)+/ \ I 

sin 012 (r)+ 4 sin fa (i+WVS-1)) * (r) J 

2m 2 r /l \ 

+> 2 x 2 4o^si r / ) ( 1 ~ 2 5 2 y5« x _ ) «r x ( ^ 12 ) f m (^X — , yx, 7 /)/< (m ’ v)v (r)=0 (4.1b) 

7, (even) ' / 


The index 7 in the above formulae represents the and in particular 
triple indices (, x, «. For example 


*/, «/ (4.2) 

A summation with respect to this index means: 


9 WWW 

2 »2 2, 2 

7 X (even) X (even) *1 (even) « x 


(4.3b) 


2*2 2 2 
V V V V 


< is a symbol which can be + or — , so that a sum 
over « consists of only these two ierms. 

(4.3a) Aside from these abbreviations, the radial equa* 
tions have been written above in a fairly explicit 
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form (with a corresponding sacrifice of elegance). 
For an arbitrary l the coupling is reasonably com- 
plicated although as we shall see in the next sec- 
tion for t small the equations simplify greatly. 
The nature of the coupling is nevertheless con- 
tained in the equations themselves via the selec- 
tion rules for the various coefficients. Regarding 
the effect of the kinetic energy coupling, a verbal 
description of the nature of that coupling is given 
in reference 7. 

The equations in the fcum ^ven~here are -ap- 
propriate to the radial coordinates r 1( r s , da with 



Also the fact that one has two functions (and 
hence two equations) for a given t and x indicates 


that the functions are in symmetric form. Hence 
the solutions may be confined to the region, say, 
ri>r* with the boundary condition 

1 =0 (4.5a) 

Jr,-r f 

and 

(4.5b) 

where. 5/&tia the-normal derivative .*n4»-<t-&nd- 
are the t^’o choices of -I- and — , being so 
selected according to (3.15) to give either the 
space symmetric (singlet) or the space antisym- 
metric (triplet) solution of the Schrodinger 
equation. 

The coupling due to the potential energy is 
given in the last term of each of the radial equa- 
tions. It can be seen that it alone is m dependent, 
and coupling coefficients l m are: 


lm(y, y\, y /) = 


4x 

(2A+1) * 


7T 2 VT J7T 

///«* 


)*£(0.*0«X £)<*•*/>'/ 


0 0 0 


sin ddddQdf 


(4.6a) 


This quadrature, which is a linear combination of integrals of the form (3.11), yields: 

T , -v _ ((l + iomfoxC^— l))(l + $0n(V2— l))(l + $,m$o«/(V2 — l))]"" 1 

lm(T,TX,rd- [«}*{<*}(<,} 

+«*«/( — 1)/+</ ) (fX — mO J f/X/)f(A —xn\h— xi) 4- </( fX — x/i | t/K,) 

+t\(C\—x—n\(f—x / ) + t(/\xj a I If— X/)] (4.6b) 


Here the «’s are to be thought as the numbers ± 1, 
and 

{*M(l+0»+f(l-0»} (4.7) 

These coefficients are subject to many restric- 
tions if they are not to vanish. Also only one of 
the C-G coefficients in the final square bracket 
will not vanish when any of these indices are 
non-zero. (When x=m = 0 all four terms con- 
tribute.) 

V. TRUNCATION AND THE zf STATE OF Hi 

The coupling in ( s and A contained in the poten- 
tial terms in principle includes a doubly infinite 
sum. If one truncates the sum in t f at 
then only a finite number of terms in the X sum 


survive. The inclusion of all X terms for a given 
L, then constitutes a given order of approxima- 
tion in the context of the partial wave theory. 

It is now our assertion that the solutions of suc- 
cessive approximation converges usefully. By 
this we mean two things; first that for a given L 
the //V)'"* 1 * 5 * from t—t, to t-L get smaller in 
some sense. Secondly in going from one L to the 
next the for a given set of indices ((, x , «) 

does not change much. 

In view of the complicated nature of these equa- 
tions this would be a very d ; fficult thing to prove 
rigorously. The heuristic argument, however, 
follows that for the one-electron diatomic mole- 
cule 1 . Namely, in the centrifugal terms, which 
are those multiplied by (r l -, -f rj"*) in the above 
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equations, there is a repulsive term proportional 
to f(f+l ) which tends to diminish the amplitude 
of each succeeding f t {r) component in the region 
of interaction. 

Since, however, there are considerably more 
terms than this f(f+l) term, the cogency of this 
argument is considerably more obscure than its 
one-electron counterpart 1 and in any event it is 


known from at least one atomic example that 
there can be resonance regions in which an addi- 
tional l component can in fact dominate its 
predecessors. 1 * 

We shal’ therefore examine one case in which 
there exist numerical results. Consider the 2,+ 
states of Hj. In zeroth order (L = C) these states 
which include the ground state, are governed by 


the equation (s ~ 1) : 


I* 


d 2 . 1 d 2 

dif l+ 7 2 drf 2 ' 


1 


• o jL\ 

sin On d0i 2 \ Sm %0i 2 / 


+ “ 2/ r i2 + 40° 


fi^ r 2 , 012)= 0 (5.1) 


This is a three-dimensional partial differential 
equation for a single function (/ss/ o <0>0)+ ). The 
expression in square brackets is manifestly in- 
variant with respect to exchange, r&zr t , and the 
singlet and triplet solutions are distinguished by 


the boundary conditions: 
[^/(r», r„ 

l/(ri, f2, 01jjr t "»r,~O 


singlet (5.2a) 
triplet (o.2b) 


Below we shall outline a systematic approximation procedure ior solving this equation. F or complete- 
ness, however, we shall give this equation with rj, r 2 , and ru as the independent variables: 


'i a 2 ,i a 2 ,2 a 2 

Sr f+n 5£“ 


ri+r?2— r 2 

r ir 2 


0 2 

dr^r-,2 


r»+r 2 2— r 2 
nr 2 


a 2 

dr 2 dr\i 



ru 


+2g°(^,r 1 )+2gc(|^a, rj ))]/(r 1 , r J( 012)= 0 (5.3) 


This form of the zeroth order equation is ap- 
propriate for solution in terms of Ilylleraas-type 
expansions. In reference 7 other ways of writing 
these equations are given, involving in particular 


asymmetric functions F 



A methodical way of solving the On equation is 
to make a relative partial wave expansion of / : 


/(fi, ri, 0«) = 2 4 , B (rir J )P,(0i2) (5.4) 

*-0 ri r s 


The Legendre polynomials P»(0, a ) are selected 
because they are the eigenfunctions of the 0« 
derivative term in (5.1). Also the interelectron 
repulsion 1/ri* has a well known expansion in 
terms of P»(0i2), so that one can easily derive an 
infinite set of two dimensional partial differential 


equation, which are essentially those that already 
have been derived in the case of the scattering of 
electrons from atomic hydrogen 11 . One can then 
truncate the sum at n = N and expect convergence 
in N in a completely analogously way as one 
expects convergence in L in the partial wave 
theory as a whole. 

Alternatively one can expand the total wave 
function in the form (considering here only the 
singlet case) : 

♦( 1 , 2 )= 2 2 2 C ni „„[R ni ^ n 1 ) 

»»»» *>»« M 

y^(fli)r^(n») + (lr=*2)l (5.5) 

If one restricts this sum to 1* and all values 
of n for each »i, then this expansion is equivalent 



GENERAL 


1345 




to (5.4). A detailed investigation of the ground 
2+ state of H 2 using (5.5) has been carried out 
by Hagstrum and Shull 14 . They have in fact 
relaxed the condition v x = -w, so that what they 
approximate in those cases is a hybrid of approxi- 
mations within the partial wpve theory. They 
have, however, examined the vi = rj=0 (spheri- 
cally averaged) approximation. This is equiva- 
lent to including one term in (5.4) within the L = 0 
(zeroth order) approximation. They obtain 
£= -2.0894 ryd. (1 ryd. = 13.6 cV) at Rab = 1.4 
Bohr radius. Tt is quite clear that the inclusion 
of all vi — Vi and with their associated m terms will 
lower this energy to below —2.1 ryd. Such a 
value is to be compared to E-— 2.34 884 ryd., 
the very accurate nonrelativistic value (including 
the nuclear repulsion) of Kolos and Roothan 1 *. 
The differences between these numbers, less than 
0.2 ryd., is very close to the difference between 
the zeroth order and exact H* f energies near their 
minimum (R^j»=2) 1 . This is just what our 
heuristic argument would lead us to expect in 
view of the fact, this being a 2 state, that the next 
angular momentum, 1=2, is two units greater 
than the included 1=0 component of the zeroth 
order problem, hi complete analogy with the case 
for the Hj + molecular ion 1 . In addition for the 
exact solutions of the partial wave approximations 
the analogy carries through regarding the upper 
boundedness of the energies of the ground and 


excited state solutions and cf their mutual 
orthogonality 1 . 
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